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We derive a master equation describing the evolution of a quantum system subjected to a sequence 
of observations. These measurements occur randomly at a given rate and can be of a very general 
form. As an example, we analyse the effects of these measurements on the evolution of a two-level 
atom driven by an electromagnetic field. For the associated quantum trajectories we find Rabi 
oscillations, Zeno-effect type behaviour and random telegraph evolution spawned by mini quantum 
jumps as we change the rates and strengths of measurement. 
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In his famous textbook, von Neumann described 'two fundamentally different types of interventions' by which a 
quantum state can change 0. These are (i) the effects of a measurement and (ii) the automatic changes associated 
with evolution via the Schrodinger equation. A subset of the latter 'intervention' is the case of open-system evolution, 
which occurs when the system of interest interacts with some other system whose properties are not necessarily of 
particular interest. Under appropriate conditions, open systems exhibit dynamical behaviour which is reminiscent of, 
and can be interpreted in terms of, the first kind of 'intervention', that is, the system evolves in a manner that reflects 
the effects of measurements being made on the system. It is the interplay between open system evolution and the 
general description of measurement performed on a quantum system that is the main focus of this paper. 

In describing the evolution of an open system, it is typically the case that we work with the reduced density operator 
of the system, obtained by tracing the joint state of the interacting systems over the degrees of freedom of the second 
system, the aim being to obtain the equation of motion of this reduced density operator, otherwise known as the 
master equation for the open system. This is done at a cost: information is lost on the correlations between the 
systems. This tracing procedure will, in general, be a highly non-trivial task except in the circumstance in which 
the open system is coupled to a bath or reservoir, often identified with the system's surroundings or environment, 
characterized as having an enormous number of degrees of freedom, and which interacts with the system over a broad 
range of energies. This interaction gives rise to irreversible dynamics associated with the loss of information from 
the system to its environment (3, |3|, U, |5| • A master equation can also be constructed if there is a fluctuating, but 
unobserved, classical element in our Hamiltonian, such as a laser phase or frequency, and it is necessary to average the 
state with respect to this random parameter ft In eithe r case, the requirements that the dynamics be Markovian, 
i.e. that the dynamics of the system be determined solely by its current state and not its previous history, and that 
the trace, Hermiticity, and complete positivity of the reduced density operator be preserved during the evolution, 
constrain the master equation to have what is known as the Lindblad form 0, 



where p is the reduced density operator of the system, H is the system Hamiltonian, bk is a system operator, and the 
7^ are positive rates. Various microscopic models of system-environment interactions yield master equations of just 
this form under the Born-Markov approximation. 

The reduced density operator in general describes a mixed state, so in a sense gives the 'average' behaviour of the 
system. But it is possible to decompose this mixed state into an ensemble of individually evolving pure states such 
that a suitable average over these states yields the density operator once again. The evolution of each such pure 
state traces out what is referred to as a quantum trajectory in the Hilbert space of the system §. The evolution 
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is stochastic in nature, i.e. each trajectory is a realisation of a stochastic process. This stochasticity takes the form 
of random, discontinuous changes in the state of the system, so-called 'quantum jumps', though under appropriate 
circumstances, these jumps will be infinitesimal and sufficiently frequent that the trajectories are continuous, and 
in general, because of this stochasticity, the quantum trajectories can be very different to each other. But what is 
truly remarkable about these trajectories is that by coupling the environment to a measuring device, it can be shown 
that the individual quantum trajectories can be understood as representing the evolution of the pure state of the 
system conditioned on the information gained about the system via the output of the measuring device. Furthermore, 
different kinds of measurements performed on the same system lead to trajectories with different stochastic properties, 
the extremes being the discontinuous and continuous possibilities mentioned above. It is also possible to generate a 
measurement record by numerically simulating the evolution of a pure state quantum trajectory; indeed the similarity 
between real measurement records and those arising from individual trajectories can be startling 

This logical pathway connecting Markovian master equations and measurement, as obtained by a microscopic 
system-environment analysis, can be inverted: instead of constructing a measurement interpretation from a master 
equation obtained microscopically, what is done here is to use measurement theory to construct as master equation 
by explicit reference to the effects of observation on the quantum system of interest. Accounting for the measurement 
record in the evolution then naturally leads us to study quantum trajectories. 

The simplest description of a measurement in quantum physics was provided by von Neumann [l|. Consider an 
observable A, represented by the Hermitian operator A, with eigenvalues {A„} and corresponding (non-degenerate) 
eigenvectors {|A n )}. A measurement of A will give one of the eigenvalues as a result, with the probability for obtaining 
the result A„ being P(X n ) = (A n |/3|A ra ). Von Neumann postulated that immediately following the measurement the 
state of a quantum system p changes to become the eigenstate of A corresponding to the measurement result : | A n ) ( A„ | . 
If the result of the measurement is not recorded then we can only state that the density operator is now diagonal in 
the basis of the eigenstates of A and that the density operator transforms as 



P 



^P(A„)|A„)(A„|. (2) 



The von Neumann measurement can be generalised to account for experimental imperfections and for measurements 
that do not correspond to a simple observable. Such generalised measurements may be described by a probability 
operator measure (POM) 0] (also known as a positive operator valued measure 0|) ^ ne elements of which, 7Tj, 
are Hermitian, positive and sum to the identity operator. Each element corresponds to a particular measurement 
outcome, and provides the probability for that particular outcome to be found by a measurement on the system. The 
probability that a measurement associated with the POM will give the result i is P(i) = Tr(/37h;). 

The POM formalism needs to be supplemented by a prescription for determining the state that the system after 
the measurement. This problem does not appear for the less general von Neumann measurement as in this case the 
system is automatically left in the measured state. The problem was addressed by Kraus who postulated forming the 
POM from effects - pairs of operators Ai and A\, the Kraus operators, such that the POM elements are 

Tt l = A\A l . (3) 

We should note that the form of the effect operators is not uniquely determined by the associated POM elements. If 
we know the form of the effect operators then observing the result i changes the density operator p to 

If the result of the measurement is not known, then the density operator will comprise a sum of such terms, each 
weighted by the probability of the associated measurement result: 

^Eno^gJpEAM!. (5) 

It is possible to further generalise by allowing each POM element to be associated with more than one effect, in which 
case 
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The effect changes the density operator for the system so that after a measurement it is left in the state 

p^J2^pA\ k . (7) 

i,k 

Note that after a particular measurement result there is no requirement for the density operator to be left in a state 
corresponding to the measurement result: TTj /Trf-ib ).pK^ 

As was mentioned earlier, there are links between measurement and master equations. Quantum trajectory reali- 
sations correspond to the system collapsing into particular states at particular times accompanied by the creation of 
some kind of measurement record by a device coupled to the environment, and thus constitute a kind of measurement- 
based interpretation of Lindblad master equations. Conversely, it is clear that interaction with a measurement device 
will cause an irreversible loss of information from the system, and thus it ought to be possible to describe continuous 
measurement by an information-loss type master equation. Attempts have been made to do this. In particular, the 
equation |3.H3| 

'p = -~[Hj}-^[d,[d,p}} = -^{H,p] + K (dpd-^d 2 p-^pd 2 ^ , (8) 

has been suggested as one which models monitoring the system observable represented by Hermitian operator O. Here 
k is a constant characterising the strength and frequency of the measurement. This equation is clearly of Lindblad 
form and expressions of this form have been derived using model couplings and environments or model measurements 
[ljl Hil Kl\ . The effect of the non-Hamiltonian term is to induce a diagonal form for p in the basis of the eigenstates 
of O. In this way superpositions of these states decohere and become statistical mixtures in much the same way as 
would result from interaction with a measurement device and subsequent tracing over this device |l8j : these states 
are commonly referred to as the pointer basis for this reason. It has also been proposed that a master equation of this 
form, with an appropriate form of pointer basis, might arise due to some intrinsic decoherence or measurement effect 
0,H(J. Master equations of this form, although widely used and much discussed 0, O^EL are n °t without their 
difficulties, however. In particular, selecting O as the position operator, corresponding to a position measurement, 
has been shown to lead to unphysical heating of the system subjected to the observation |2^ . 

It is not clear, however, how we can connect master equations of this form with the description of the measurement 
process presented earlier. Nor is it clear that such a comparison is possible in all cases. This question is important if we 
wish to be sure that the evolution does indeed model a measurement and, more importantly, to quantify the accuracy 
of the measurements and follow the acquisition of information in individual evolutions. Clarifying the connection 
between Markovian master equations and measurements is the main objective of this paper. 



DERIVATION OF THE MASTER EQUATION 

In this section we provide a simple but comprehensive derivation of a master equation describing the continuous 
monitoring of a quantum system. Our starting point is to use the Kraus operators to describe the effect of the 
measurements on the system state. We assume that the continuous monitoring takes the form of a sequence of 
measurements which take place instantaneously and randomly in time, but at an average rate R. Consider a short 
time interval At. During this time the probability that a single measurement will occur is RAt and we suppose 
that the time interval is short enough so that the possibility that two or measurements occur can be neglected. If a 
measurement does occur then the density operator transforms as 

p{t)^p{t + At)=Y,U[t)^. (9) 

i 

If no measurement occurs then the density operator follows the normal Schrodinger evolution and if At is sufficiently 
small then we may evaluate this change to lowest order in At: 

p(t) -> p(t + At) = p(t) - ~[H, p(t)]At. (10) 

The average evolution results from adding these two, weighted by their probability of occurence. To lowest order in 
At we find 



p(t + At) = (1 - RAt)p{t) - ^[H, p(t)]At + RAt^2Aip(t)Al 

i 



(ii) 
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If we take the limit At — > then we obtain the measurement master equation 



-[H,p] + R 



(12) 



This equation, which is clearly of Lindblad form, is the main result of our paper. A similar approach has also 
been suggested by Stenholm and Suominen to form a master equation describing a continuously monitored system 
[2^ |. More generally, we may need to include measurements with POM elements of the form JJjJl. In such cases the 
measurement master equation will take the form 



h 



[H, p] + R 



ik 



(13) 



We have applied this equation to describe quantum friction in terms of effective measurements of a particle's position 
and momentum associated with collisions with molecules forming the surrounding medium [24L |25| . It is also possible 
to derive from this equation the more general master equations of Lindblad form [H] and ^ which will be published 
elsewhere. 

The measurement master equation describes evolution due to both weak and strong measurements: the parameters 
governing measurement strength are hidden within the Kraus operators. The measurements may also be frequent 
or infrequent and it is possible to play these parameters off against one another. We might expect a high rate of 
weak measurements to give similar evolution to a lower rate of stronger measurements. This does indeed appear to 
be the case for the evolution of the density operator and we will find in our discussion for a two-level atom that the 
measurement master equation depends only on a single rate formed from the observation rate R and the a parameter 
describing the strength of the measurement. When we examine the quantum trajectories, however, the observation 
rate and the measurement strength have quite distinct effects on the evolution. Indeed, quantum trajectories can 
exhibit wildly different types of evolution. We should note that the quantum trajectories have a greater significance 
than the master equation as we will have access to the measurement record. Individual trajectories enable us to 
simulate both properties of the system and the associated measurement record. 



MEASUREMENT MASTER EQUATION FOR A DRIVEN TWO-LEVEL ATOM 

The simplest quantum system has just two distinct states and this suffices to illustrate the construction and 
application of our measurement master equation. We consider two energy levels of an atom |1) and |2), resonantly 
driven by an electromagnetic field 0, @- An alternative physical realisation is a spin-1/2 nucleus subjected to a 
constant magnetic field oriented in the z-direction and a radio-frequency field with its magnetic field pointing in the 
x-direction 26] . It is convenient to describe the two-level atom in terms of the three Pauli operators 

o x = |1)<2| + |2)(1| 
a 2 = i(|l)(2|-|2)(l|) 

a 3 = |2)(2|-|1)(1|. (14) 

For the driven two-level atom these correspond respectively to the real and imaginary parts of the atomic dipole and 
to the atomic inversion. In a suitable interaction picture, the Hamiltonian describing the resonant interaction between 
the atom and the driving laser field (which we describe classically) has the form 

H = -—cn. (15) 

This interaction causes the probability for the atom to be found in the state |1) to oscillate in time with the Rabi 
frequency O The corresponding fully quantum problem has the classical laser field replaced by a single quantised 

field mode and leads to much more complicated dynamics |27j . 

We suppose that during the Rabi evolution the single two-level atom is interupted by a series of instantaneous 
imperfect generalised measurements of the energy level. These are represented by the POM elements 

7r 1 =p|2)(2| + (l-p)|l)<l| 

7r 2 =p|l)<l| + (l-p)|2)(2|, (16) 
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associated with the measurement suggesting that the atom was in state 1 or state 2 respectively. Here < p < 1/2 
is the probability that there is an error, i.e. the probability that a measurement performed on the state |2) gives 
the result 1 or vice versa. When p = the measurements are perfect, and determine, without error, the state of the 
system, when it is small we have strong measurements and when p is close to 1 /2 the measurements are weak. At 
p = 1/2 both POM elements are proportional to the identity and the measurements tell us nothing about the state 
of the system. 

In order to find an appropriate measurement master equation we need the effect operators. These are not uniquely 
determined by the form of the POM but, for the sake of simplicity, we can choose the following as the most simple 
forms consistent with the POM elements: 



A 1 = Vp|2)(2| + ^r^|l)(l|=Al 

A 2 = ^\l)(l\ + ^T~p'\2)(2\=Al (17) 

A measurement of this kind can be realised in nuclear magnetic resonance by means of coupling the spin under 
consideration to one of its neighbours j^- These can be substituted directly into the general measurement master 
equation to give 

= y[oi,p] + — -p-^pj (a 3 pa 3 -p). (18) 

It is clear that the master equation depends on the rate and measurement error probability only through the combina- 
tion 7 = -f (V 1 — p — y/p) . A similar dependence of an effective rate, 7, on the rate of observation and the accurac 
of the measurements can be found in models in which the system under observation is coupled to a meter system 
Furthermore, this equation is equivalent to eq. (JBJ describing a measurement of 03, as <rf = I. We will find, however, 
that varying R and p leads to dramatically different quantum trajectories induced by the measurements. We can 
solve equation (|18|) by writing the density operator as 

p = - (/ + u{t)a x + v(t)a 2 + w(t)a 3 ^j . (19) 

This leads to coupled equations for the components of the Bloch vector (u,v,w) — ((<ti), (<t 2 ), (03)) S@ : 

u = —2ju 

v = Clw — 2jv 

w = -Qv. (20) 

Solution of these is straightforward and gives 

u{t) = u{Q)e- 2 ^ 

v(t) = w(0)e" 7 * (cos(n'i) - ^ 7 sin(fi't)) + w(0)e-^ ^ sm(n't) 

w(t) = w{0)e~^ (cos(fi't) + sin(Q't)) + v^e^^- sin(ft't), (21) 

where fl' = (f2 2 — 7 2 ) 1 / 2 is the Rabi frequency as reduced by the action of the measurements. This evolution 
is reminiscent of that which occurs when the laser driving the atom is subject to phase jumps |29| or frequency 
fluctuations [sfl ]. 

In fig. H we plot the evolution of the inversion w(t) as a function of the dimcnsionless time ttt. We see that the 
Rabi oscillations are damped out and that the inversion tends to a steady state value of zero. This is a reflection 
of the fact that the steady state for the atom is the fully mixed state p(oo) = corresponding to zero values for 
u, v and w. We also see that that the period of the Rabi oscillations is increased by the measurements through the 
presence of 7 in f2'. If 7 > f2 then we will not have any Rabi oscillations and the evolution of the inversion is then 
reminiscent of over-damped harmonic motion. 
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FIG. 1: A plot of atomic inversion as a function of time for 7 = 0.1414J! 



MEASUREMENT-RESULT DEPENDENT EVOLUTION 



The most obvious benefit of our approach is that we can easily follow the evolution of the atom conditioned on the 
results of the measurements. We can do this by implementing the change Q after a measurement, but it is simpler 
and more natural to employ the quantum trajectories method [8|, l3ll |32( to simulate individual realisations of an 
experiment by calculating the evolution of the state vector In this method, we break time into short intervals 
of length St and then use a random number to determine whether to simulate a measurement (which we perform 
with probability R5t) or otherwise to evolve according to the Schrodinger equation with the Hamiltonian 115|) . If 
a measurement is performed then the result, 1 or 2, is chosen by means of a random number in accord with the 



probabilities (A\Ai) and (A r 2 A2) respectively. If the measurement result is 1 then the state changes as 

AM) 



with a corresponding change for the measurement result 2. 



U/2 



(22) 



Weak measurements 



If the error probability p is close to 1/2 then the quality of our measurements will be poor and the effect operators 
will be close to the identity operator: 



A 2 



y/l-p+y/P) 1 + [Vi-P-Vp jO-3 
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(23) 



The unit operator part leaves the state of the system unchanged, and only the small part proportional to 03 alters the 
state. Figure |2] shows a typical trajectory of the state of the system in the regime where the measurements are very 
weak. The value of the overall decay rate 7 is the same as in fig. Q which shows decaying Rabi oscillations, but in fig. 
El the Rabi oscillations never decay. The measurements do almost nothing to the state, and so the state just evolves 
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FIG. 2: Atomic inversion as a function of time for weak measurements. Parameters are p — 0.49, R = 20fi, so the equivalent 
decay rate 7 = 0. 14140 

like an undamped two-level atom. The principal change from undamped Rabi oscillations is that the Rabi period is 
increased slightly. We can understand this in terms of the effect of the measurement on the state. After a measurement, 
the action of one of the operators A\ and A2 is to reduce the real and imaginary parts of the dipole u and v while 
the other increases it. The measurement result inducing a reduction in the dipole is always the more likely outcome, 
however, so that on average a measurement reduces the dipole by the factor 2y/p(l — p) w 1 — 4(1 — 2p) 2 . This leads 
the sequence of measurements to induce the dipole to decay at the rate 2j and, following eq. I|2U|) . the reduction in 
the size of the dipole is responsible for a deceleration of the evolution of the inversion w. As the measurements are 
frequent, the gross features of the evolution of the inversion are determined by this average behaviour. The decay of 
the oscillations which appears in the solution of the measurement master equation (fig. ^) is a consequence of the 
randomness of both the results of the measurements and of the time at which they occur. This causes a dephasing 
between the Rabi oscillations in different realisations. 

In fig. |3 the measurement is of slightly better quality and we see that the Rabi oscillations are more strongly 
perturbed by the measurements. It is clear that some of the measurements make a significant discontinuous change 
in w. As the measurements become stronger (corresponding to ever smaller values of p) these discontinuous changes 
will eventually become quantum jumps between the states |1) and |2). In fig. 0]we see the effects of more frequent 
but weaker measurements: the discontinuities occur much more frequently but the changes associated with each are 
correspondingly smaller. In the measurement master equation there is a trade-off between increasing R and reducing p 
and it is possible to do both whilst leaving 7 and the evolution of p unchanged. In the individual quantum trajectories, 
however, increasing R and decreasing p leads to qualitatively different evolution. It is also apparent that the Rabi 
period is increased beyond its value in fig. [2] This is a consequence of the fact that the effects associated with the 
measurement results decrease the values of u and v by a greater amount and hence slow further the Rabi oscillations 
of w. 

Perfect measurement and the Zeno effect 

If the probability for error is small (p ~ 0) then we are in the regime of strong measurements. If the rate of 
measurement is low then we will find regular Rabi oscillations interrupted by changes of state either to |1) or |2), 
followed by the resumption of oscillations. This is reminiscent of the trajectories found in the study of resonance 
fluorescence in the presence of spontaneous emission If the measurements are very good (p = 0) and are carried 
out very frequently then we find a quantum Zeno- type behaviour [4| . The measurement repeatedly collapses the system 
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FIG. 3: Atomic inversion as a function of time for weak measurements. Parameters are p — 0.36, R = 1.4140, so the equivalent 
decay rate 7 = 0. 14140 




FIG. 4: Atomic inversion as a function of time for weak measurements. Parameters are p — 0.49, R = 258.80, so the equivalent 
decay rate 7 = 18.300 

back onto one of the energy eigenstates before it has time to evolve away significantly. Rabi evolution is seemingly 
inhibited completely and the system remains in the excited state for a long time. The system occasionally jumps from 
one state to the other after an average time which depends on the measurement rate and the Rabi frequency. This 
is more likely to occur immediately after there has been a relatively long time between successive measurements, as 
the state then has had time to evolve significantly away from one of the eigenstates. A measurement which finds the 
system in the other eigenstate then has a reasonable chance of occurring. 
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For perfect measurements the POM elements and Kraus operators are 



7T1 

7To 



|2)(2| 



Ao. 



(24) 



Let us suppose that the state evolves under the influence of the Hamiltonian 1|I5|) for a short time 6t between two 
measurements. If the state after the first measurement was |2) then in this time it will evolve to 



= VT^\2)+e\l) 



(25) 



where e « flSt/2 <C 1. The probability that the measurement of this state finds the system in state |1) is (tti) = e 2 , 
which is very small. The most likely result is that the measurement will find the atom in the state |2) and this will be 
the state of the atom after the measurement. In this way the measurements suppress the Rabi evolution leading to 
stability of the states |1) and |2). This frozen motion is punctuated by occasional jumps between the states leading 
to the random telegraph behaviour illustrated in fig. Very occasionally the interval between the measurements 
will be sufficient for the coherent evolution to cause the inversion to depart noticeably from the stable values of ±1. 
These are just visible in fig. [S] The departures of the inversion from ±1 and the jumps themselves are suppressed as 
the measurement rate increases. This is because the probability that the measurement will induce a jump depends 
on the square of the time elapsed since the preceding measurement. In the limit of infinitely frequent measurements 
the state then remains in either the initial excited or ground state indefinitely. 




-0.5 
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FIG. 5: Atomic inversion as a function of time for perfect measurements. Parameters are p = 0, R = 100S7, so the equivalent 
decay rate 7 = 50£1 



Imperfect measurement, small jumps and the weak Zeno effect 



If the measurement error probability is small but finite we get an entirely different type of behaviour. The Rabi 
oscillations are still suppressed, but the system undergoes a short-period random telegraph. Furthermore the inversion 
also exhibits small and short-lived excursions from I or -1 (fig. 0). These filaments are associated with small 
measurement-induced quantum jumps, not from the excited state to the ground state, but between intermediate 
superpositions of these states, beginning and ending near the maximum absolute values of the inversion. The existence 
of such filaments is the signature of what we refer to as the weak Zeno effect: a measurement can be associated with 
beginning a transition but the next measurement can stop it. What appear as full jumps from I to -1 and vice versa 
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FIG. 6: Atomic inversion as a function of time for strong measurements. Parameters are p — 0.16, R — 70.86f2, so the equivalent 
decay rate 7 = 18.30f2 




FIG. 7: Detailed evolution near a change in inversion in fig. (JSJ 



are in reality the effect of several of these mini-jumps occurring quickly, as is borne out by the expanded-timescale 
versions of this figure (figs. 0and|HJ|. 

Consider the effect on a general state = a|l) + /3\2) of sequences of measurements. We assume that, as in 
the figure, the rate of monitoring is sufficiently rapid that evolution does not change the state significantly between 
successive measurements. After one measurement the system is in one of the states \ip\) oc j4i or \ip2) oc ^l^): 

lw = V(i- P )ilp+ P i/y (^ |1> + ^ |2> ) m 
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FIG. 8: Detailed evolution near a filament in fig. © 



I '02 



■p\a\ 



(27) 



After two measurements there are four possible states, given by the four possible measurement sequences AiAj with 
»,i = l,2: 



1 



11/ 



((l-p)a\l)+pl3\2)) 



\, •,.,} = |V2i) = ; = = (jpjl=p~a\l) + v^^^/3|2>) = |V->) 



1 



y/(l-p)\/3\ 2 +p\t 



(pa|l) + (l-p)/3|2)) ; 



(28) 
(29) 
(30) 



where the square of the normalising denominator of each state is the probability for the associated sequence of 
measurement results to occur. The two states corresponding to different successive measurements occurring are the 
same as the original state. Different successive measurement results have no nett effect: the inversion can be induced 
away from +1 or —1 and then jump back by the same amount following the next measurement. The probability of 
this occurring is low (~ p) for strong measurements but such jumps are apparent in fig. [H| An expanded view of one 
such sequence is given in fig. |H1 One of the other two possibilities is much more likely, however, with either both 
measurements results being 1 or 2. 

If the state is an almost equal superposition \a\ ~ |/3| then this will evolve rapidly towards one of the energy 
eigenstates. The probability of three or more identical successive measurements is again significantly enhanced. The 
overall effect is that approximately equally weighted superpositions of |1) and \2) are unstable. 

A qualitative understanding of the dynamics can be obtained, therefore, by considering the effect of sequences of 
measurements on states of the form of eq. (12511 which has a — e ~ and f3 = \[\ — e 2 ~ 1. After the first measurement 
the state will become 



1^0 
1^2) 



l-e 2 2 



el 



VPl 2 > 



v^|i) 



p)e 2 +p(l-e 2 ) 



yr^yr^|2) 



v / J5 e 2 + (l-p)(l-e 2 ) 



|2>. 



(31) 
(32) 
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Thus measurement result 2 leads to no significant change to the system, but measurement result 1, which occurs with 
the low probability e 2 + p, increases the probability that the system is in state |1). This amounts to a small jump 
from state |2) downwards, and is the beginning of one of the filaments in fig. The amount by which the system 
jumps is determined by the relative sizes of e and yfp. 

The time between jumps is small, so again we can disregard the evolution, which amounts to the assumption that 
e 2 -C P- We therefore only need to consider the effect of a second measurement on \ipi) 



, = (1-P)e\l) +P VT^\2) e|l) +P |2) 
^(l-p) 2 e 2 +p 2 (l-e 2 ) ^e 2 + p 1 

= ^^ £|1) + ^^^ |2) * |2). (34) 
^ P (l-p)e2+p(l-p)(l-e2) 

Thus if the second measurement result is 2 the system returns to the initial state in line with the general state result. 
If the second measurement is 1 then the proportion of state |1) is further increased. The ratio of the probabilities 
of the measurement outcomes 1 and 2 is (e 2 +p 2 )/p, an increase on the ratio for the first measurement (which was 
~ e 2 +p). Thus if a filament does start, then its continuation is more likely than was its start. 

The analysis in the previous paragraph shows that an occasional mini-jump down associated with the unlikely 
measurement result means that for the next measurement the jump downwards is more likely, and can be of the same 
order of likelihood as the jump back up. If the jump is back upwards then the result is one of the mini filaments seen 
in fig. Any combination of jumps down and up which ends back at the excited state results in a filament. The 
greater the number of excess number of down-jumps the longer the filament. In some cases the filament extends below 
an inversion of zero so that the state of the system contains a greater proportion of the ground state than excited. 
The situation is now reversed and jumps down become more likely. When this occurs and the filament reaches the 
ground state (inversion -1) what appears to be a full jump has occurred. These jumps seem to be single jumps straight 
down when the system is viewed on a coarse timescale, but finer timescales reveal that they consist of several of these 
mini-jumps, both down and up. The rate at which these full jumps occur is significantly higher than for the pure Zeno 
jump rate, as the probability of several consecutive mini-jumps is high compared with the probability of Zeno-type 
single jumps for perfect measurements. 

When the system gets close to the ground state the size of the mini-jumps decreases until they are comparable in 
size to the typical Rabi evolution between jumps. Then a tiny amount of Rabi evolution introduces a small amplitude 
associated with state |2) and the whole mini-jump process begins again, but with the roles of the excited and ground 
state reversed. The overall behaviour of the system over long times is then a short-period random telegraph, in which 
the jumps from the excited to the ground state consist of several mini-jumps. 

The mini-jumps persist for quite large measurement error probability p and for significant Rabi evolution between 
measurements. As p decreases, however, the initial jump size increases, until for p — the jumps are directly from 
close to the excited state to the ground state, and the Zeno situation is reached. With decreasing p the filaments and 
the weak Zeno effect disappear as the likelihood of several mini-jumps occuring becomes small. 

Also of note is that some of the realisations pictured here are performed with parameters p and R such that the 
overall decay rate is the same, for example, figs. 0]andElwhose overall decay rate is 18.3 times the Rabi frequency. In 
fig. 01the measurement error probability is 0.49, and the measurement rate is 258.8 times the Rabi frequency, whereas 
for fig. the paramenters are 0.16 and 70.8. Thus wildly different behaviours can occur for the same decay rate. 
Another less dramatic example of different evolutions for the same decay rate is given by figs. [5] and [3] for which the 
damped evolution associated with the solutionof the master equation is depicted in fig. ^ These examples show that 
for a master equation to describe continuous monitoring well it must have a two-component decay rate. 



CONCLUSIONS 



In this paper we have derived a Lindblad master equation which describes the evolution of a system under continuous 
monitoring. Instantaneous measurements of the system are assumed to occur at an average rate. The equation is 
based on Kraus's effect description of the measurements, which provides information about the state of the system 
after the measurement process. The equation describes evolution under both weak measurements, where the state (no 
matter what it is) is hardly altered by an individual measurement, and strong measurements, which can change the 
state significantly. In contrast to previous equations derived for this purpose, whose only free parameter is the decay 
rate, the measurement master equation decay rate contains independent measurement rate and strength parameters, 
so quantum trajectory analyses based on this equation will be qualitatively different even though the overall decay 
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rate may be the same. We will further develop the links between the measurement master equation and general 
Lindblad master equations in future work. 

We have applied our measurement master equation to the resonantly-driven two-level atom, deriving the master 
equation based on monitoring the energy eigenstate of the atom. We have performed a quantum trajectory analysis 
of the system, and find qualitatively different types of behaviour for different rates. For weak measurements the 
individual quantum trajectories resemble the undamped Rabi oscillations. The principal effect of the measurements 
is to increase the Rabi period and to randomise the phase of the oscillations. For perfect measurements we can obtain 
the suppressed evolution of the well-known Zeno effect. The system remains in an energy eigenstate, |1) or |2), for a 
long time, after which it jumps instantaneously to the other eigenstate. This behaviour continues, with the system 
inversion performing a long-period random telegraph type evolution. 

For strong measurements, where the probability of measurement error is small but finite, the system again performs 
a random telegraph evolution. The 'flat' parts of the random telegraph signal are filamented: there are very short 
departures from the eigenstates which are induced by the measurements. These are due to smaller mini quantum 
jumps, both up and down, between intermediate superpositions of the two states. The filaments result from departures 
from the eigenstates which end up where they started. The full jumps in the telegraph signal are seeded by, and consist 
of mini-jumps, again both up and down, between intermediate superpositions of the two states, but ending in the 
other eigenstate. The mini quantum jumps, or weak Zeno effect, are an unusual phenomenon, particularly, as we are 
not monitoring the states into which the system jumps. Under these cicumstances they can only occur for imperfect 
measurements of the state. It is difficult to see how they could be found systematically using other approaches in 
which the measurement strength and rate are not separate parameters within a single decay rate. By starting with 
the measurements, and the associated Kraus operators, the independence of these two parameters is clear from the 
outset in our approach. 
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